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We present a simple and physically compelling boundary condition regularization scheme in the 
framework of effective field theory as applied to nucleon-nucleon interaction. It is free of off-shell 
ambiguities and ultraviolet divergences and provides finite results at any step of the calculation. Low 
energy constants and their non-perturbative evolution can directly be obtained from experimental 
threshold parameters in a completely unique and model independent way when the long range 
explicit pion effects are removed. This allows to compute scattering phase shifts which are, by 
construction consistent with effective range expansion to a given order in the CM momentum and 
are free from finite cut-off artifacts. We illustrate how the method works in the 1 5o channel for the 
One Pion Exchange potential. 
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Effective field theories (EFT) have been successfully investigated in recent years in the context of hadronic and 
nuclear physics. Their main ingredient has to do with the occurrence of scale separation between long and short 
distance physics, making the development of a systematic power counting possible. Since the original proposal of 
Weinberg's [ij to make a power counting in the potential many works have followed implementing such a counting J2, 
H, 0, w ith finite cut-offs or proposing a counting in the renormalized S-matrix 0, Q which has also been pursued [J. 
Both Weinberg and Kaplan-Savage- Wise schemes can be understood as perturbative expansions about infrared fixed 
points @ (see also Ref. [13). In any case, converge improves under certain conditions pd| . According to Ref. |l2| a 
hybrid counting involving also the chiral limit should be invoked (see also Ref. ^H). For a recent and more complete 
review on these and related issues see e.g. Ref. 0] and references therein. 

Much theoretical insight has been gained by analysing how short and long distance physics separate for the One 
Pion Exchange (OPE) interaction in the singlet 1 Sq channel where the scattering length, uq — — 23.7fm, is much 
larger than the size of the potential l/m^ = 1.4fm. The renormalization of NN interaction in this channel has been 
studied several times in the literature. In Ref. an elegant subtraction method has been developed to construct a 
finite T matrix for contact, i.e. zero range, interactions added to OPE. Renormalization is indeed achieved by taking 
the subtraction scale to be much larger than any other mass scale and checking for independence of results in this 
limit. The resulting description of the 1 Sq is only valid to very low energies, requiring for inclusion of derivative terms. 
Unfortunately, the method cannot be easily extended in that case. These derivative interactions can be included within 
a cut-off regularization [13 , Il7j and dimensional regularization 0] . In this latter case a three-parameter fit can be 
achieved with no explicit two pion exchange contribution. For the pionless theory, though, the inconsistency between 
both regularization methods after renormalization has been pointed out [l9|. Momentum space treatments based 
on the Lippmann-Schwinger equation appear more natural from a diagrammatic point of view within a Lagrangian 
framework and allow explicit consideration of nonlocal potentials. In practice, however, the long range potentials 
used in NN scattering are local, and for those the analysis of renormalization in coordinate space may be simpler. In 
addition, the Schrodinger equation is a second order operator and boundary conditions define a complete solution of 
the problem in the whole space both inside and outside the boundary. This is equivalent to a sharp separation between 
the interior and exterior region. This property is naturally formulated in coordinate space for a local potential. 

Although the idea of using boundary conditions for NN scattering is a rather old one (see e.g. Ref. 1201 and references 
therein), there have been recent works in this regard motivated by the developments within EFT |2lLl22tl23| . Actually, 
it has been shown that in the absence of long range forces a low momentum expansion of the potential within 
EFT framework for the Lippmann-Schwinger equation is completely equivalent to an effective range expansion and 
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also to an energy expansion of a generic boundary condition at the origin in coordinate space for the Schrodinger 
equation. If a long range OPE potential is added we will show below that due to the short distance Coulomb nature 
of this potential the origin must be reached continuously from above R — > , R > (i.e. excluding the point R = 0), 
in harmony with known theorems on self-adjoint extensions of Schrodinger operators |24| . 

In this paper we analyze precisely how the energy dependent boundary condition must change as we move the 
boundary radius for fixed energy to achieve independence of physical observables such as scattering phase shifts. By 
doing so we are effectively changing the Hilbert space since the wave function in the outer region is defined only from 
the boundary to infinity. An advantage of this procedure is that we never need to invoke off-shellness explicitly; at any 
step we are dealing with an on-shell problem. In addition, we work directly with finite quantities and no divergences 
appear at any step of the calculation when the boundary radius is taken to zero from above. 



II. VARIABLE PHASE EQUATION WITH BOUNDARY CONDITIONS 

The reduced Schrodinger equation for including OPE in the channel for NN-scattering with CM momentum k 
reads 

-u'l{r) + U{r)u k {r) = k 2 u k (r) , (1) 
together with the asymptotic condition at infinity 

Uk(r) -> sin(fcr + S(k)) . (2) 

The OPE potential in the 1 So channel reads 

U{r) = -i^L e - . (3) 

167r/^ r 

Where Af/v is the nucleon mass, the pion mass, the pion weak decay constant and qa the nucleon axial 
coupling constant. In the numerical calculations below we take Mn — 938.92 MeV, = 93MeV , = 138 MeV 
and gA — 1.25. Our unknowledge of the interaction below a certain distance scale R is parameterized in terms of a 
boundary condition at the matching point r = R, 

u' k (R)-L(k,R)u k (R)=0. (4) 

In general, this boundary condition depends both on the boundary radius R and the momentum k. The value of R 
separates the whole space into two disjoint regions, an outer region where we assume the interaction to be given by 
OPE potential, and an inner region where interaction is regarded as unknown. 

The boundary condition at R 7 Eq. Q) has a simple physical interpretation. If we switch off the long range piece 
U(r) above the scale R, then the phase shift due to the short distance physics below the scale R is given by 

= L(k,R) = kcot(kR + 6(k,R)). (5) 

u k (R) 

It is interesting to see what kind of equation satisfies the short distance phase shift, 6(k, R), as we steadily move the 
boundary radius R for a fixed momentum k. Using Schrodinger 's equation at the boundary r — R we get the variable 
phase equation, 

d5{k R) = _ 1 gin2 + 
dR k 

The obvious condition at infinity must be satisfied 

lim S{k, R) = S(k) (7) 

R— >oc 

Thus, Eq. © describes the evolution of the phase shift as we go down to lower distances, assuming that both the 
long distance potential and the physical phase shift are known. Regardless of whether or not the potential we are 
considering is realistic at very short distances 1 one can extrapolate the long distance potential to the origin and 



1 Two Pion Exchange becomes comparable to OPE at about the distance of r = 1.5 fm. So, any extrapolation of Eq. J^J with OPE below 
1.5 fm should not be considered realistic. 
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define the zero range OPE-extrapolated phase shift 

6 s (k) = lim S(k,R). (8) 

The precise manner how this limit is built depends on the OPE potential, Eq. J3J|, and will be analyzed below. Eq. © 
is well known in potential scattering ( for a review see e.g. Ref. |25|1. but it has always been used assuming the trivial 
initial condition Ss(k) = lim/j„>o S(k, R) = 0. 



III. LOW ENERGY EXPANSION OF THE BOUNDARY CONDITION 

The former variable phase equation, Eq. @ can be cast in a more convenient form by defining the variable 
K —matrix, 

K{k,R) = k cot S(k,R), (9) 

yielding 



dK(k,R) rr/ „s \ T ,,, ^smkR , 1 

V J — - = U(R) K(k,R) — hcosfci? 

dR k 



(10) 



At low energies, however, it can be conveniently parameterized as an effective range expansion, which carries over to 
the variable phase 

kcotS(k,R) = \— + \r Q {R)k 2 + v 2 {R)k i ■ ■ ■ (11) 



a (R) 2 



one has 



U(R) (ao - Rf (12) 



dctQ 
~dR 

4rr , n J 1 /r R 1 \ n ( ^ R \ ( 1 r f v 2 1 R , 1 



These equations have to be supplemented with some initial conditions ao(i?o)j r a(Ro) an d ^2(^0) at a given boundary 
radius, i?o- If we take the initial boundary radius, Ro = the set of equations, |Q , (| 1 31) and l|14|l express the evolution 
of the low energy parameters at short-distances when the long distance potential is switched on. Conversely, if the 
initial boundary radius is taken to infinity they offer a possibility to determine the short-distance low energy parameters 
from the total ones by downwards evolution in the variable R when the long distance potential is switched off. Notice 
the very appealing and natural hierarchy in the previous equations; while the distance evolution of the scattering 
length ao is autonomous, the remaining low energ y pa rameters ro, v 2 , etc. depend on the previous ones. The relation 
to the renormalization group approach of Ref. [9LTTo| will be discussed elsewhere |23 ■ 

Before presenting the numerical results l|12|). I|13l) and (|14f> we analyze first the short and long distance behaviour. 
At short distances R << l/m n the OPE potential behaves like the Coulomb potential. Eq. (|12fl can be easily solved 
in two cases, ao << R and ao >> R. While in the first case we get 

a (R) = a (R ) - &™*% N (R 2 - R 2 ) , a « R (15) 

in the second case one solution behaves as 

r m _ Qo(^o) 16tt/^ 1 

a ° [ ] 1 + ao (R o )^0^\o g (R/R o )" 9 2 A mlM N \o g (R/R ) ' °° <K " ^ 

where R < R a << l/m v . As we see, ao(-R) goes to zero very slowly and with a' (R) — > —00 at short distances, 
which in momentum space corresponds to the ultraviolet limit. The resemblance of Eq. (|16fl with asymptotic freedom 
is striking. It is easy to see that the first case, Eq. i|15|) . corresponds to selecting the regular solution at the origin, 
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FIG. 1: Evolution of the scattering len 1 So NN-threshold parameters ato(R) (in im),ro(R) (in fm) and V2(R) (in fm 3 ) from 
the asymptotic values at infinity (which we take in practice Roc = 10 fm ) when OPE effects are removed down to the origin. 
oto = — 23.73 fm and ro = 2.68 fm and V2 = — 0.48fm 3 . Solutions of Eqs. 1)12113114^ are labelled as "exact". The extrapolated 
values at the origin when OPE effects are removed are as,o = 0, rs.o = 4.08 fm , and vs,2 = 0.43 fm 3 . We also show some 
approximations for cea(R). OPE means One-Pion- Exchange only and corresponds to integrate Eq. (1 1 Lit from the origin to 
infinity with the boundary condition Qo(0) = 0. SDE means short distance expansion as given by Eq. I16H . LDE correspond 
to a long distance expansion, Eq. I17H (natural case) and Eq. 1181 (unnatural case), respectively. 



whereas Eq. (|16|l is the generic case, which always contains an admixture of the irregular solution. Obviously, the 
regular case is exceptional and for that particular situation one can integrate from the origin starting with the trivial 
initial condition <5(fc,0) = up to infinity. The result corresponds to a pure OPE interaction, with no short-distance 
interactions. The important thing to note here is that no matter what the initial value of arj was at infinity (except for 
the exceptional case discussed before) , removing one-pion exchange goes into the same value at the origin, as implied 
by Eq. (|16[l ). This also means that any small deviation of the ao(Ro) at small distances results in huge variations at 
infinity. Thus, removing OPE results in a extreme fine tuning of the low energy parameters at short distances. 

We analyze now the long distance behavior. Clearly, when R >> l/m n we have ct' (R) = 0, Eq. I|12|). and we 
approach quickly the asymptotic value ao(oo). For such long distances we can always use perturbation theory to 
solve the equations backwards. For scattering lengths which are natural, i.e. ao ~ I/t^ we may neglect ao(R) with 
respect to R and get 

/>oo 

a (R)-a = U(R)R 2 dR+ ... (17) 

Jr 

For unnatural scattering lengths, cvo >> l/m T we make the opposite approximation, and get 

1 1 



a (R) a J R 



U(R)dR+... (18) 



The previous equations (|17fl and 118fl hold irrespectively of the strength of the potential, provided R is sufficiently 
large. Similar approximations for the remaining low energy parameters will be discussed elsewhere [27J. 

The numerical evolution of as t o(R) and r$(R) starting with the experimental values, a — —23.739 fm, r = 2.68 fm 



5 



and V2 = —0.48 fm 3 down to the origin according to Eqs. 1|12I13I14|) is shown in Fig. Q) 2 . We also show the perturbative 
estimate in the case of large and small scattering lengths based on a long distance expansion Eq. I)17|l (natural case) 
and Eq. (|18|1 (unnatural case), respectively, as well as our short distance estimate, Eq. I|16() . In the case of cto(R) 
we observe a huge change from infinity down to the origin, although remains unnatural, ao(R) » R. Numerically 
we confirm our theoretical expectation that 0^0(0) = (see Eq. I|16|) ). This simply means that the renormalized 
contact interaction becomes arbitrarily small as the OPE potential is switched off. This is, however, not the case for 
the renormalized derivative interaction, as expected from our estimate, Eq. (|16[1 . Our numerical values extrapolated 
to the origin are 

as,o = ao(0 + ) = (19) 
r Sfi = r (0+) = 4.08fm (20) 
v s .2 = v 2 (0 + ) = 0.43 fm 3 (21) 

This is the initial condition which, in principle, has to be supplemented in Eqs. H12I13I14|) in order to get the 
experimental results (see also discussion below). The work of Ref. [2(| uses a two Yukawa model to extract the 
short-distance low energy parameters. This is done by fitting the data and then switching off the OPE contribution, 
yielding 015.0 = — 1.72 fm, r$fi = 1-60 fm and vs.2 = — 0.024fm 3 . In Ref. [IJ an attempt to determine the short- 
distance parameters based on the three Yukawa model yields as,o = — 3.38 fm, rs,v = 2.60 fm and vs,2 = 0.313 fm 3 . 
The short-distance scales in that calculation are R p = 2/m p = 0.46 fm and R a = 2/m a = 0.80 fm. For that range 
we get a,? n = —3.6, — 2.21fm, rs i0 = 2.7, 3.1 fm and v$,2 — 0.18, 0.26 fm 3 respectively, in qualitative agreement with 
Refs. [l2,l29- Note, however that our way of determining the short-distance low energy parameters does not require 
any specific model at short distances. 



IV. ^o-PHASE SHIFT 



Once the short distance parameters are known one may compute the phase shifts to any order of the approximation 
in a fc 2 expansion of the initial condition without any additional parameter fitting by integrating Eq. I|10|) upwards 
with a suitable initial condition at a short distance radius, 

K s (k) = k cot S s (k) = — + \r^ s k 2 + v 2 , s k A + ■■■ (22) 

The standard way of proceeding is to determine the low energy constants or equivalently the short distance parameters 
directly from a fit to the data and then recompute the threshold parameters. A further advantage of avoiding a fit 
is that one can prevent spurious and/or multiple minima; our solution is essentially unique. Morevover, since by 
construction at a given order in the k 2 expansion the low energy behavior of the phase shift is reproduced up to the 
same order in k 2 , the possibility of getting even slightly different threshold parameters due to a fit in the intermediate 
energy region is precluded. Actually, our procedure would coincide with the standard one, if the fit was carried out 
in the region where an effective range expansion holds (k < 60MeV if v 2 is included). 

Due to the fact that the origin is a fixed point for the running scattering length, i.e. ao(R) — > for R ^ regardless 
of the value of ao(o°), Eq- EH, one must integrate the equations from very small distances upwards, using the value 
of ckq(R) at that distance. It is important to realize that a tiny mismatch in the value of oq close to the origin results 
in a complete different value of ao and also of the phase shift at infinity 

In Fig. @ we show the results for the phase shift depending on the number of terms kept in the low energy 
expansion at short distances. Our results exhibit a good convergence rate. For comparison we also depict the effective 
range expansion results without explicit pions, which is expected to work at low energies only, and corresponds to 
make R$ — > 00 in our approach. As we see, the effect of introducing pions always improves the results. This can 
be fully appreciated at NNLO, where ER does a poor job above CM momenta ~ lOOMeV, but explicit OPE effects 
enlarge the energy range up to about ~ 140MeV ~ m n . where we expect explicit two pion exchange contributions to 
start playing a role. 



2 In practice results are insensitive for long distance cut-off of Roo = lOfm. In the case of the short distance cut-off we can go down to 
Rs = 0.0001 fm without much effort but results are fairly insensitive to the short distance radius already at ijg = 0.1 fm, where we 
have a§ q = — 0.9867fm, rog = 3.819fm and «js = 0.397fm 3 . For shorter distances Eq. 1161 provides an accurate estimate for ao(R). 
Taking larger values of R B builds in finite cut-off effects. Actually Rs >> l/m^ corresponds exactly to effective range expansion. 
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FIG. 2: Predicted phase shifts according to Eq. 1101 when OPE potential is switched on and the initial condition is a low energy 
expansion of the K-matrix at short distances, (See Eq. I22H in the main text). LO means keeping as,o only, NLO keeping 
as,o and ro,s and NNLO keeping as,o, ro,s and V2,s- The short range parameters are directly determined by evolving the low 
energy parameters from their experimental values cto = — 23.73 fm and ro = 2.68 fm and V2 = — 0.48 fm 3 . ER-LO, ER-NLO and 
ER-NNLO corresponds to a pure effective range expansion keeping cto only, as,o and ro, cto, ro and i>2 respectively. OPE-only 
corresponds to OPE without short-distance contributions. No further fit is involved. Data are the PWA from Ref. |28| . 



An interesting point to note at this stage is that if as,o = with other short-distance low energy parameters 
fixed, we would inevitably get 6s(k) = mr, as deduced for instance from Eq. If we solve the variable phase 

equation with that condition at R — up to R — R^ >> l/m T we get the result (also shown in Fig. © for 
comparison) corresponding to a regular OPE with the regular boundary condition Uk(0) — instead of the mixed 
boundary condition of Eq. (0J at R = 0. The puzzle is resolved by realizing that the limiting procedure in the 
boundary condition and the solution do not commute; the limit R — > + implies 5'(k,R) — > oo whereas starting at 
R = requires 5(k,R) ~ R 2 producing instead a bound derivative 5'(k,R) ~ R (see Eq. ©). This discontinuous 
dependence of the boundary condition on the boundary radius at R = agrees with rigurous theorems on self-adjoint 
extensions of Schrodinger operators (see e. g. Appendix D of Ref. [24|). 



V. CONCLUSIONS 



In the present paper we have analyzed the renormalization of the OPE interaction in the presence of contact and 
derivative interactions of any order for NN scattering. In order to do that we have derived an equation for the 
evolution of an energy dependent boundary condition in coordinate space as a function of the boundary radius. The 
resulting equation shares many properties with renormalization group equations and can be interpreted in terms of 
the phase shift produced by eliminating OPE from infinity to the boundary radius, which eventually is taken to 
zero. Two advantages can be deduced from this framework: no divergences appear and there is no need to consider 
off-shell extrapolations. This allows to set up equations for the running low energy parameters as a function of 
the boundary radius. Using the experimental values for the low energy parameters, which correspond to an infinity 
boundary radius, we extract in a unique and model independent way the corresponding short-distance parameters. 
Our numerical values agree with other determinations based on specific models for the short-distance interaction. 
As we get closer to the origin we find a fixed point structure, triggered by the non-vanishing contribution of the 
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irregular solution. This requires a fine tuning of the short-distance low energy parameters. After that we integrate 
the running phase shift upwards and determine without any additional fit the 1 Sq phase shift. The OPE plus contact 
and derivative interactions to NNLO is able to describe the 1 5o phase shift up to CM. momentum of about 140 MeV, 
which coincides with opening of the two pion exchange left cut. Above that momentum explicit two pion exchange 
effects should set in. 

The results presented in this paper are very encouraging and suggest several improvements and extensions. Explicit 
Two Pion Exchange contributions are expected to contribute significantly at about 1.5 — 2fm, so our results should 
not be considered realistic below that scale, or equivalently above CM momenta of about 100 — 150MeV, as it seems 
to be the case. In addition, our description should be enlarged to include all partial waves. Work along these lines 
will be presented elsewhere p7| . 
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